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Abstract 

Classical systems of coupled harmonic oscillators are studied using the Carati- 
Galgani model. We investigate the consequences for Einstein's conjecture by con- 
sidering that the exchanges of energy, in molecular collisions, follows the Levy type 
statistics. We develop a generalization of Planck's distribution admitting that there 
are analogous relations in the equilibrium quantum statistical mechanics of the 
relations found using the nonequilibrium classical statistical mechanics approach. 
The generalization of Planck's law based on the nonextensive statistical mechanics 
formalism is compatible with the our analysis. 

Key words: 

Planck's law, nonequilibrium, superstatistics 



1 Introduction 



It was observed that for blackbody radiation [I] the mean energy U of sta- 
tionary electromagnetic waves is a function of frequency v of these waves: at 
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temperature T we can write 



U(u ->■ 0) ->■ fc B T and U(u ->■ c») ->■ 0, 



(1) 



where fee is the Boltzmann constant. This contradicts the energy equipartition 
law of classical statistical mechanics in that U is frequency independent. The 
important contribution of Planck came when he discovered that one could ob- 
tain the experimental result of Eq. (1) by treating the energy of the stationary 
electromagnetic waves as a discrete spectra. Planck's law opened the way for 
quantum mechanics. 

In such a way, Einstein [2] gave an important contribution by noticing that in 
the canonical ensemble the fluctuation of energy of the blackbody radiation is 
compatible with Planck's law when it is, in terms of its variance <te, expressed 
in terms of the mean energy U through the relations 



where (3 = 1/(/cbT), e = hv (h is Planck's constant) and N is the number of 
quantum particles. Einstein imagined the first relation as being exactly a type 
of general thermodynamic relation, while the second must have a dynamical 
character, being, in principle, deducible from a microscopic dynamics. 

Recently, Carati and Galgani [3ll4] have investigated classical systems of cou- 
pled harmonic oscillators with the procedure of Jeans-Landau- Teller dynamic 
of molecular collisions [5ll6] . The principal results of Carati-Galgani are that: 
i) relaxation times to equilibrium are nonuniform and depend on the internal 
degrees of freedom of the considered system; ii) Situations of nonequilibrium 
statistical mechanics very far from equilibrium are described, in a first order 
approximation, by the analogous relations found in quantum equilibrium sta- 
tistical mechanics; iii) as conjectured by Einstein, Eq. (J3J) is a consequence of 
the dynamics. 

It is important to remark that Planck's law considers a real equilibrium state 
of the quantum statistical mechanics, while within the framework of Carati- 
Galgani, a similar relation to Planck's law was obtained using nonequilibrium 
classical statistical mechanics. 




(2) 



and 




(3) 
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The purpose of the present work is to study classical systems of coupled har- 
monic oscillators along the line of Carati-Galgani. We studied three different 
aspects of this framework. First, we analyzed the consequence for Einstein's 
conjecture by assuming that the exchanges of energy, in molecular collisions, 
follows the Levy type statistic for high enough frequencies. The Levy statistics 
for the energy exchanges was recently obtained by Carati et al. [7] , considering 
the dynamics over a finite time in the Jeans-Landau- Teller model. Secondly, 
admitting that there are analogous relations in the equilibrium quantum sta- 
tistical mechanics of the relations found here using the nonequilibrium classical 
statistical mechanics approach, we will explore some possible consequences in 
quantum radiation theory. Third, we show that the generalization of Planck's 
law [8] within the framework of the recently introduced superstatistic formal- 
ism [9fT0] is compatible with the our analysis of Einstein's conjecture. 




The organization of this paper is as follows. The Carati-Galgani model is 
briefly explained in Sec. II. Our generalization is presented in Sec. III. Discus- 
sions about the consequences of our results are presented in Sec. IV. Conclud- 
ing remarks are given in Sec. V. 



2 The Carati-Galgani model 

A good level of interest has been focused on the problem of estimate energy 
exchanges between vibrational and translational degrees of freedom in the 
dynamic of molecular collisions. One of the simplest and main model used to 
study this problem was published by Carati and Galgani [3]H] . It is based on 
the Jeans-Landau- Teller approach [5]. A possible example of a physical system 
that reproduce this formalism follows. 

Let us imagine a one-dimensional system involving two particles P and Q in 
a line. P is attracted by a linear spring to the fixed origin through a harmonic 
oscillator potential. Q interacts with P through an analytical potential. In a 
first approximation, the variation of energy 6e in a simple collision between 
P and Q for the oscillator P of frequency v is given by the addition of a drift 
term and others dependent on the fluctuation of the oscillator's initial phase. 
It is easy to show [5] that after t collisions an oscillator has energy 



eo is the oscillator's initial energy that is considered to be very small in order 
that the motion of P and Q be decoupled, ipj is the phase of the oscillator in 
the jth collision and r\ is evaluated as a function of z/, which is known to be 





(4) 
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decreasing exponentially with v (for an analysis of rj see [6]). 

Let us consider a assembly of oscillators, averaging over the phases (fj, in order 
that 

t 

(et) = (e )+t77 2 + 277 ^( v ^-Icos(^_i)). (5) 

3=1 



where (...) is the phase average. Considering the phases to be independent and 
uniformly distributed, for large t, the average of energy and the variance are, 
respectively, given by 

(e t ) = (e )+t V 2 (6) 



and 

a 2 et = 2eobH) + (/i 2 > . (7) 



Here a\ and (n t ) correspond to, respectively, the variance and mean of the 
exchanged energy (e t — to). Observe that the functional form of the variance 
is independent of t and rj. 

Assuming N independent identical oscillators and considering the total energy 
(E t ) = J2j(l^t), we immediately obtain for the variance a\ = Na 2 and mean 
U = N(n) of the exchanged energy 

(a E ) 2 = 2,uU + ^ (8) 



where <; is the initial action per oscillator [3] . This result represents essentially 
the application of the central limit theorem, that arise from a large number 
of independent contributions of the oscillators [TT] . 

Equation (jSJ) is analogous to Einstein's relation obtained in the Bose-Einstein 
quantum statistical mechanics theory. Despite the fact that this relation has 
been obtained considering the classical statistical mechanics very far from 
equilibrium, the approach recovers Einstein's conjecture, i.e., proving it to be 
a consequence of the dynamics. 
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3 The generalization of the Einstein's conjecture 



Recently, Carati et al. [7J showed that there are systems which present qualita- 
tive differences between the statistics obtained using the equilibrium distribu- 
tion and that obtained through the dynamics over a finite time. Considering 
the dynamics over a finite time of the Jeans-Landau- Teller model, they ob- 
tained that the exchanges of energy follow the Levy type statistics for high 
enough frequencies. An immediate consequence of this result, based on the 
Levy flight random walk [T5], is that (fi^) is infinite. 

We notice that in the Carati-Galgani model it is assumed that the oscillators 
have low energies and the collisions are separated from each other so that the 
oscillator phases are completely independent. Here, we must be able to answer 
what happens if the (fif) is not finite? In this way, we can no longer write Eq. 
([7]) and, as consequence, Eq. ([8]) is not valid either. 

It is well known that for (fj,^) — > oo the theoretical basis is the Lindeberg-Levy 
central limit theorem [13]. Since we are considering the case of low energies 
for the oscillators, using the same steps as were used to obtain Eq. flTJ, it is 
straightforwardly verified that 

al = 2e (jj) + e i2 - a) (^ Q ) 1 < a < 2 , (9) 



where e is a constant that has dimension of energy and is suitably chosen by 
dimensional requirements. 

Finally, for N independent identical oscillators, we find 



a % = 2wU + eW^- I (10) 



which generalizes the dynamic proposal of Einstein and will be our hypothesis 
to get the generalizations of Einstein's conjecture and Planck's law. 



4 Possible relations in equilibrium statistical mechanics 



Let us now explore the possibility of there existing analogous relations in the 
equilibrium quantum statistical mechanics of the relations found here using 
the nonequilibrium classical statistical mechanics approach. 

We start the process using Eq. ([2]) in conformity to Einstein's vision. In this 
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case, using Eq. ( !T0|) in Eq. fl2]) we have 



It is straightforward to find the energy distribution 

^-iV{ £(2 _ a)[e( ^_ i] } ^, (12) 



in order that we can consistently rewrite Eq. (Tl2"j) as 
Ev 



We have introduced U = — H = 2q and d = with 1 < d < oo. 
We assumed £ oc Hv that leads to the expression e = zHv with z a pure 
dimensionless number. For d = 1, H coincides with Planck's constant h, and 
the U associated with Bose-Einstein quantum statistic is reobtained. 

In the same way as the Planck blackbody radiation has its origin in the Bose- 
Einstein statistics, here obtained using d — 1, we can take other values of 
d, so that we can obtain new physical systems such that Eq. ffl3|) may be 
applied. By following along standard lines, after straightforward calculations, 
the density of energy per unit volume is obtained as 

Vl d (u) = g(u)U = (14) 



z' 



l c 3[ e /3Hv/d _ iy 



where we use the standard density of states of the photons (g{v) = 87r^ 2 /c 3 ). 
Fig. [T] presents fid(z/) versus energy Hv for typical values of d and (3 = 1. Note 
that if d increases then Vt^v) increases. Although ^(0) has a finite value for 
v — > 0, it is null for d < 3 while it is infinity for d > 3 (see inset). The Qd is 
modified significantly when d varies. The fid are illustrated in Figs. [2] and [3] 
for, respectively, d = 2 and d = 4 and various values of f3. For d fix, we note 
small changes for fid when varies. 

In the Hvfi -C 1 regime 



8vri/ 3 - £/ 
[zHY^c 3 



toM - 7ZT7UZ^(dK B T) d (15) 



which generalizes the Ray leigh- Jeans law. The generalization of Wien's law is 
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Fig. 1. The density of energy per unit volume versus energy for d = 1,2,3,4 and 
(3 = 1. (Inset shows a closeup of part of the graph for v — > 0) 




Fig. 2. The density of energy per unit volume versus energy for /3 = 0.9, 1.0, 1.1 and 
d = 2. 

obtained in the Hv(3 ^> 1 region 

~ "m e KbT > ( 16 ) 



and the Stefan-Boltzmann law for the total emitted power per unit surface is 
attained from the usual one, given by j = <7dT 4 , where the d-dependence is 
present only in the prefactor cr^- 
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Fig. 3. The density of energy per unit volume versus energy for /3 = 0.9, 1.0, 1.1 and 
d = 4. 

5 Souza-Tsallis generalization from the Einstein's conjecture 



In fact, one of the characteristics of Eq. (|T3|) is that it reproduces the results of 
Souza and Tsallis [8]. Observe that for U = Hun we can obtain the particle- 
number distribution 



n 



[ e PHv/d _ iy ' 



(17) 



Adopting the formal theory of statistical mechanics, it is straightforward to 
find the variance of an ensemble of n particles as being 



o\ = n + n a 



which generalizes the bosonic Landau's relation |15j . 

Ref. [8] obtained the Eqs. (JT7J) and f|T8l) heuristically based on the nonextensive 
statistical mechanics [16] and Beck-Cohen superstatistics [9fT0] . They were 
based on the discussion of the differential equation dy/dx = —aiy — (a q — ai) y a 
(with y(0) = 1), whose particular case a = 2 corresponds to Bose-Einstein 
statistics. Our formalism describes the Souza-Tsallis generalization from the 
Einstein's dynamical fluctuation viewpoint. 

Finally, we see that Eq. (fT7|) can be written as 

oo 

n = Y,9(hd)e-? E % (19) 

i=0 
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where we obtain the energy spectrum 



Ei(d) 



His 



(i + d), 



(20) 



d 



and its degeneracy 



g(i,d) 



r(i + d) 



(21) 



r (d) r (i + 1) ' 



T(x) being the gamma function. The spectrum is made of equidistant levels, 
like that of the quantum harmonic oscillator, with a constant energy difference 
between successive levels equal to Hv/d. 



6 Conclusions 

Summarizing, we introduced a generalization of Einstein's conjecture con- 
sidering an approach along the line of the Carati-Galgani model of molecu- 
lar collisions. We obtained Planck's distribution compatible with the Beck- 
Cohen superstatistics. Furthermore, establishing a metaequilibrium very far 
from equilibrium for the system studied, and admitting that there are analo- 
gous relations in the equilibrium quantum mechanics, we generalized Planck's 
law. 

For instance, the effects in the quantum radiation theory were obtained the- 
oretically, as a first order approximation, of a classical system very far from 
equilibrium [4]. We emphasize that we are working with classical models of 
molecular collisions and not quantum models for blackbody radiation. Under- 
standing the relationship between these two points of view is an open problem. 
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